Energy protection arguments fail in the interaction picture 
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Spin Hamiltonians with degenerate ground states are one potential system for the storage of 
quantum information at low temperatures. Trapped ions can be used to simulate the dynamics 
of these Hamiltonians, but the coherence-preserving properties will be lost. This illustrates that a 
quantum simulation performed in an interaction frame will not thermalize with its environment. 



-j^-, I. INTRODUCTION 

o: 

1 The preservation of quantum information is difficult. The system must be protected from noise due to its environ- 
ment and from errors induced by control and measurement. In comparison, classical information is robust. Usually, 
, there is a large energetic/kinetic barrier separating the two states relative to the temperature. A common example is 
I ^ ' the use of magnetic domains for computer memory 

Motivated by the idea of magnetic domains and topologicalparticles, a number of ground state degenerate quantum 
' systems have been proposed as error free memories [1, UI3, ESQ- ^^^h case, the logical qubit is the degenerate 
ground state of many physical qubits in the presence of a Hamiltonian, H . The environmental noise interacts with 
I— I the qubits independently and H is constructed such that all single qubit errors raise the energy. Consequently at 
I a temperature lower than the energy gap, multi-order processes are required to decohere the system. Additionally, 
although the lowest order processes do not preserve the exact state, they often preserve the superposition over large 
subspaces. This is equivalent to quantum error correction where the distance between encoded states allows one to 
detect and correct up to n errors. The ground states of these Hamiltonians can be used to derive quantum codes for 
^ ■ quantum error correction 0,0. However, the goal is to avoid active error correction and let the energetics maintain 
Q^l the coherence. In these systems, the error rate is suppressed as l/A'^ where A is the energy gap between the ground 

■ state and k depends on the size of the system and the details of the Hamiltonian. 

' The Hamiltonians used for protecting coherences are traditionally associated with condensed matter physics but have 
^ . now been proposed for atomic and molecular physics using arrays of ions, atoms, and molecules [lOl [Til IT^ . Atomic 
ions have already been shown to have long lived qubits relative to the frequencies of these effective Hamiltonians fl3| . 
f — ■ Is it possible that these methods can be used to extend these lifetimes even longer? 

' A problem arises because the effective Hamiltonians are in a rotating frame relative to the Bohr frequency of the 
' physical qubit. Although the energy gap may protect the system against noise of frequency lower than A, it will not 

■ suppress the natural spontaneous emission. The utility of such a scheme then depends strongly on the kinetics. 

' In this paper, the effects of spontaneous emission in these systems is examined by looking at a model two qubit 
C , system. We find that the interaction picture is equivalent to a system in an environment with negative energy modes. 
We then examine the implications for energy-preserving Hamiltonians, quantum simulations, and adiabatic quantum 
computations in the interaction picture. 



Oh' 



o 



X 



^ . II. 2 QUBITS COUPLED BY XX + YY 

A. Model System 

The model systems is described by the two qubit Hamiltonian 



H^v = ^{W2 + YiY2) (1) 

where X and Y are Pauli spin operators. The eigenstates of the Hamiltonian are — "^(|01) " IV-'2) — 

[00) , IV'a) = 111) , and IV'4) = 75 (|01) + |10)) with eigenenergies Ei = -hJ, E2 = E3 = 0, and E4 = hJ. 

Although there is no protected qubit since the ground state is not degenerate, has the property that any single 
qubit Pauli error on jV'i} raises the energy by at least hJ. Therefore, we expect that the state 1-01) will be preserved for 
temperatures lower than hJ/kt. We can quantify this intuition by assuming each qubit interacts with an independent 
bath of harmonic oscillators and calculating the transition rate from \ipi) to all other states. 
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The Hamiltonian of the system and the environment as 



H = Hxy + ^ huJka\f.aik + ?i ^ aikzZi{a\f^ + Qik) + aikxXi{a\f^ + Qik) + aikyYi{a\^. + aik) (2) 



i,k 



where aiki are the interaction strengths of mode k with qubit i along the axis j. 

We can transform to an interaction picture with respect to the bath Hamiltonian, Ho = H^, = ^ fiijjka\^aik. The 
Hamiltonian in the interaction picture in general is H — exp{iHot/h){H — Hq) exp{~iHot/fi). For this case, we find 

H = Hxy + hJ2a^kzZ,iale'^''' + a.ke-"^"') + a^kxX,{a\^e'^''' + a^fee"*"''*) + a.kyY.ial^e'^"' + a^fee"*"''*) (3) 

i^k 

Fermi's golden rule is used to calculate the first order transition rate from the ground state to the other states 



assuming a generic spectral density p{uj) and thermal occupation, n{uj,T) 



e^bT - 1 



We find that 



r = ATrazi2 J f p{2J)ni2J,T) + AiriaxiJf + ay{jf)p{J)n{J,T) (4) 

where a(w) is the interaction strength averaged over the modes of frequency lu. This is the heart of these protection 
schemes. Because we are in the ground state, energy must be absorbed. At low temperatures, a linear decrease in the 
temperature will lead to an exponential decrease in rt(J, T) and the transition rate. At zero temperature n(J, T)=0 
and the rate of excitation from the ground state vanishes to first order. In this case, because the ground state is 
non-degenerate, the higher-order rates will also be zero. 



B. Effective Hamiltonian with Trapped Ions 



An effective XX or YY coupling between two ions has been demonstrated by virtually exciting the shared motional 
mode [3, llBl- Although, the same motional mode cannot be used to generate an XX + YY coupling, it is possible 
to use the center of mass mode of the two ions for the XX coupling and the stretch mode for the YY coupling. This 
yields again the Hamiltonian, H^y The coupling strength J/2 = rf'Vl^ /6 where ry is the Lamb-Dicke parameter, Q, is 
the Rabi frequency of the qubit transition and 5 is the detuning of the laser from the relevant motional sideband. J 
is limited to be smaller than the ion trap motional frequency uJs- In practice, uJs/'i'K is typically a few MHz. 

Implicit in the above description is that H is in the rotating frame or interaction picture defined by the Bohr 
frequency of the physical qubit, v. For concreteness, consider the optical qubit of ^°Ca"'" where the state is |0) 
and the iS'1/2 state is |1) pJi]. For this quadrupole transition, v/2-k = 411 THz and 7/27r = 0.16 Hz. 

The state will not be preserved. Spontaneous emission will drive the system towards \(f)2) = |11). Although, 
\(f)i) is the ground state of H^y, it is an excited state of the system Hamiltonian Hs = Figured] We start with 

the same system-bath interaction as above but we then transform to an interaction picture defined by the uncoupled 
system and bath Hamiltonians Hq = Hs + Hb- The transformed Hamiltonian is then 



H^Hxy + hY, {a^k.Z,ale'^''' + [a^k. + ia,ky)a-^e-'•'\ale'^''* + ake-'^"*) + H.c}j (5) 

i.k 

where cr_i = {Xi + iYi)/2 is the raising operator on the zth qubit and H.c. is the Hermitian conjugate. 

Notice that in the interaction picture, the bath frequencies have been effectively shifted by a frequency ±1/, Figure 
[21 This is equivalent to the bath having modes of negative energy. In this situation, spontaneous emission can lead 
to heating of the system. We again calculate the transition rates 

r = AT:az{2jf p{2J)n{2J,T) + 2T:{ax{v + Jf + ay{v + Jf)p{v + J)n{v + J,T) 
+ 2T:{ax{v - Jf + ay{v - Jf )p{v - J){n{v -J,T) + 1) 

(6) 

At zero temperature, we have a non-zero transition rate 

F = 2TT{ax{iy-Jf + ay{v^jf)p{v~J). (7) 
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Interaction Picture Schrodinger Picture 



FIG. 1: Comparison of the effective system spectra in the interaction picture and Schrodinger picture. The ground state of 
the effective Hamiltonian is not the ground state of the total system, |3). The figure is not to scale. For optical qubits, 
v/ J « 10*, and for hyperfine qubits, z^/J ~ 10'^. 

and as expected F « 7 when J is small compared to as is the case in ion traps. We see that spontaneous emission 
drives the system out of the ground state of the effective Hamiltonian and the state has a lifetime equivalent to the 
excited state of a single ion. 

Theoretically, we can examine the case where the ions emit into the same photon mode. In this case, photon 
emission from the state is prevented by destructive interference, ^sub = 0. This is perfect Dicke subradiance 
[l7| and an example of a 1-dimensional decoherence-free subspace [3 13 • this situation, the applied Hamiltonian 
would reduce the errors associated with low frequency errors. In practice, ions are a few /^m apart from each other 
and subradiance is a weak effect, ^sub = 0.997 [20] ■ 

For optical ion qubits, the applied Hamiltonian does not reduce the errors. Quantum information ion trap experi- 
ments often use hyperfine states instead of optical qubits [2ll | . A hyperfine qubit is a naturally energetically protected 
subspace with respect to electronic transitions. Spontaneous emission via a magnetic transition is possible but is neg- 
ligible on the time scale of an experiment. The main source of environmental decoherence in an ion trap experiment is 
uncontrolled magnetic fields (2^ . This decoherence can be minimized by encoding into clock states |l3l . [2^. Applying 
the Raman beams necessary to generate Hj.y introduces a new error pathway due to Raman scattering [2j|. Again, 
this error is not suppressed by the applied H^y because the energy shift associated with J is small compared to the 
frequency of the scattered photons. 

III. IMPLICATIONS OF SPONTANEOUS EMISSION FOR THE INTERACTION PICTURE 

A. No energy protection 



A long-range quantum compass model has been proposed to protect qubits in an ion trap [25| . Contrary to Ref . [25| , 
spontaneous emission is not suppressed even though the emission occurs at a slow enough rate such that eigenstates 



system bath 
a) natural 



system bath 
b) effective 



FIG. 2: Comparison of the system and bath spectra. We assume a physical bath of positive energy modes, a) The expected 
spectrum for two qubits naturally coupled by H^y. If the temperature is less than the energy gap to the excited state, excitations 
will be exponentially suppressed, b) In an interaction frame, where the effective interaction is Hxy the system spectrum is the 
same. However, the bath modes are shifted to effective negative energy. Consequently, a negative energy Boson can be emitted 
allowing for the spontaneous excitation of the system. 

of the interaction Haniiltonian will decay to eigenstates of the interaction Hamiltonian. Here we examine the effects 
of spontaneous emission for the 4 qubit system. In this case, the system is equivalent to the compass model [26| and 
the Hamiltonian is 



The Hamiltonian commutes with X1X2 and X3X4. The eigenvalues of these operators label four four-dimensional 
subspaces, Aij. The logical |0) and |1) are the ground states of the subspaces A++ and respectively. It could 
be possible that spontaneous emission, although exciting the ground state, may preserve the superpositions between 
the subspaces. Unfortunately, spontaneous emission connects the subspace A±± to both subspaces and A |_. 

The eigenstates, li^j), for all the subspaces are computed and the rates between the states are calculated to be 
proportional to {i^ — {Ej — Ei)/h)^\qji\'^ where {qji) is the sum of the individual quadrapole moments evaluated for 
the eigenstates i and j. We note that {Ej — Ei)/h is on the order of J and in ion traps always much less that 1/ (even 
for hyperfine qubits). Therefore, in our calculations we approximate {v — {Ej — Ei)/h)^ w v^. 

The population dynamics are calculated starting from the logical state |0). The population randomly walks among 
the eigenstates with a rate of the natural spontaneous emission. The state that does not spontaneously emit, |1111), 
is not an eigenstate of the effective Hamiltonian; the resulting steady state is a dynamic equilibrium with equal 
population in all states. Figure [3] shows the result of this calculation. The two measures used are the population in 
the state |0) and the extent of the population in the subspace A++. We see that in a time 5/7 (5 s for **°Ca+) the 
populations are almost completely mixed. Compared to the spontaneous emission rate from the |0) state of a single 
qubit there is no improvement. 



H = nj{XiX3 + X2X4 + Y1Y2 + Y3Y4) 



(8) 
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FIG. 3: (Color online) The compass model in the interaction picture does not suppress errors due to spontaneous emission. 
The decay of the population of the logical |0) (solid line) is comparable to the decay of the excited state in the physical qubit 
(dotted line). The same is true if we consider the population of the subspace ^4++ associated with the logical |0). The applied 
compass model Hamiltonian drives the system from the bare ground state (all ions in the IS) state) and the result to zeroth 
order in J/i/ is a completely mixed state. 



B. Quantum simulations improperly thermalize 

A quantum simulation uses a controlled quantum system to model another quantum system 12711 . The principle of 
these devices have been demonstrated in the interaction picture using nuclear magnetic resonance |28l. [29l . [30 | . Trapped 
ions are seen as a promising candidate to simulate various magnetic systems 10] including the compass model described 
above. Unlike state preservation, a quantum simulation could benefit from noise. Classical techniques for describing 
noisy quantum systems are often computationally more expensive. For example, calculating the dynamics of a system 
under the Lindblad equation (3l| requires quadratically more resources when compared to the Schrodinger equation. 
Although in the limit of large noise a quantum system can be efficiently classically simulated (32i] , there should be a 
range of noisy quantum simulations that are still more efficient than a classical simulation. 

In the interaction picture the effective noise model is often not physical. In Section fill Al we calculated that the 
dynamic equilibrium of the compass model with spontaneous emission was the completely mixed state. Although the 
physical bath is at zero temperature, the effective temperature of the system is infinite. For the two qubit example, 
the final state of the system is non-Boltzmann distributed with respect to the effective Hamiltonian. Ion trap and 
nuclear magnetic resonance quantum simulations require that the simulation be over before the state decoheres. 

Quantum simulations often employ an initialization by adiabatic evolution [lol.[30|. Adiabatic quantum computation 
[3^ can be viewed as a special case of a quantum simulation where the ground state of the simulated Hamiltonian 
solves a computational problem. The ability of the system to follow the ground state is determined by the speed of 
the evolution and the minimal gap between the ground and first excited state. Although a noisy environment can ruin 
the adiabatic algorithm when the gap becomes small |3 4|] , t he thermalization of the system may erase these errors by 
recooling to the ground state as the gap becomes large [35| . Unfortunately, in the interaction picture a thermal state 
is not achieved. Instead, spontaneous emission drives the simulation / adiabatic quantum computation to a dynamic 
equilibrium that may not be the desired state. 
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IV. CONCLUSIONS 



The simulation of a system with Hamiltonian, He, is the mapping of the dynamics of onto another system. 
Experience with how the inital system interacts with a physical environment leads us to assume that the simulated 
dynamics will decohere in a similar way. This intuition is incorrect. Often the simulated dynamics are in a rotating 
frame or interaction picture. In this frame, the bath modes appear to have negative energy and the ground state of 
He is not a steady state even at zero temperature. 

Trapped atomic ions have been used to demonstrate the error-reducing techniques of quantum error correction 
[3^ and decoherence-free subspacesfs^l- However, the technique of constructing a multi-body Hamiltonian with a 
naturally preserved ground state is incompatible with how multi-qubit gates are implemented in ion traps. As a 
result, the effective Hamiltonians do not yield the coherence properties associated with the natural Hamiltonians 
proposed in 0, S i, i, S, 0, 111 . The realization of coherence-perserving Hamiltonians will require the construction 
of sophisticated condensed matter quantum information systems [ssj . 



V. ACKNOWLEDGMENTS 



K.R.B. thanks Dave Bacon for many useful discussions. This work was supported by the Georgia Institute of 
Technology. 



[1] C. A. Ross, H. I. Smith, T. Savas, M. Schattenburg, M. Farhoud, M. Hwang, M. Walsh, M. C. Abraham, and R. J. Ram, 
J. Vac. Sci. Technol. B 17, 3168 (1999). 

[2] A.Yu. Kitaev, Annals of Physics 303, 2 (2003). 

[3] J. P. Barnes and W. S. Warren, Phys. Rev. Lett. 85, 856 (2000). 

[4] D. Bacon, K. R. Brown, and K. B. Whaley, Phys. Rev. Lett. 87, 247902 (2001). 

[5] C. Wang, J. Harrington, and J. Preskill, Annals of Physics 303, 31 (2003). 

[6] L. B. loffe, M. V. Feigel'man, D. Ivanov, M. Troyer, and G. Blatter, Nature 415, 503 (2002). 

[7] D. Bacon, Phys. Rev. A 73, 012340 (2006). 

[8] D. Gottesman, Phys. Rev. A 57, 127 (1997). 

[9] P. Aliferis and A. W. Cross (2007), eprint quant- ph/0703230 
[10] D. Porras and J. L Cirac, Phys. Rev. Lett. 92, 207901 (2004). 
[11] A. MicheU, G. K. Brennen, and P. ZoUer, Nature Physics 2, 341 (2006). 

[12] A. M. Rey, L. Jiang, M.Fleischhauer, E. Demler, and M. D. Lukin (2007), eprint 'cond-mat/0703108', 

[13] C. Langer, R. Ozeri, J. D. Jost, J. Chiaverini, B. DeMarco, A. Ben-Kish, R. B. Blakestad, J. Britton, D. B. Hume, W. M. 

Itano, et al., Phys. Rev. Lett. 95, 060502 (2005). 
[14] K. M0lmer and A. S0rensen, Physical Review Letters 82, 1835 (1998). 

[15] P. C. Haljan, K.-A. Brickman, L. Deslauriers, P. J. Lee, and C. Monroe, Phys. Rev. Lett. 94, 153602 (2005). 

[16] F. Schmidt-Kaler, S. Guide, M. Reibe, T. Deuschle, A. Kreuter, G. Lancaster, C. Becher, J. Eschner, H. HafTner, and 

R. Blatt, J. Phys. B: At. Mol. Opt. Phys. 36, 623 (2003). 
[17] R. Dicke, Phys. Rev. 93, 99 (1954). 

[18] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 (1997). 

[19] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. Lett. 81, 2594 (1998). 

[20] R. G. DeVoe and R. G. Brewer, Phys. Rev. Lett. 76, 2049 (1996). 

[21] D. J. Wineland, C. Monroe, W. M. Itano, D. Leibfried, B. E. King, and D. M. Meekhof, J. of Res. of the National Inst, of 
Standards and Technology 103, 259 (1998), http://nvl.nist.gOv/pub/nistpubs/jres/103/3/cntl03-3.htm. 

[22] D. J. Berkeland, J. D. Miller, J. C. Bergquist, W. M. Itano, and D. J. Wineland, J. Appl. Phys. 83, 5025 (1998). 

[23] P. C. Haljan, P. J. Lee, K.-A. Brickman, M. Acton, L. Deslauriers, and C. Monroe, Phys. Rev. A 72, 062316 (2005). 

[24] R. Ozeri, C. Langer, J. D. Jost, B. DeMarco, A. Ben-Kish, R. B. Blakestad, J. Britton, J. Chiaverini, W. M. Itano, D. B. 
Hume, et al., Phys. Rev. Lett. 95, 030403 (2005). 

[25] P. Milman, W. Maineult, S. Guibal, L. Guidoni, B. Dougot, L. lofle, and T. Courdreau Phys. Rev. Lett. 99, 020503 (2007). 

[26] J. Dorier, F. Becca, and F. Mila, Phys. Rev. B 72, 024448 (2005). 

[27] S. Lloyd, Science 273, 1073 (1996). 

[28] D. Cory, R. Laflamme, E. Knill, L. Viola, T. Havel, N. Boulant, G. Boutis, E. Fortunato, S. Lloyd, R. Martinez, et al., 
Fortschr. Phys. 48, 875 (2000). 

[29] C. Negrevergne, R. Somma, G. Ortiz, E. Knill, and R. Laflamme, Phys. Rev. A 71, 032344 (2004), arxiv:quant-ph/0410106. 
[30] K. R. Brown, R. J. Clark, and I. L. Chuang, Phys. Rev. Lett. 97, 050504 (2006). 
[31] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 



7 



[32] D. Aharonov and M. Ben-Or, in Proceedings of 37th Annual Symposium on Foundations of Computer Science (FOCS) 
(IEEE Computer Society Press, 1996), p. 46. 



[33] E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser (2000), e-print |quant-ph /OOP 1 106 
[34] M. S. Sarandy and D. A. Lidar, Phys. Rev. Lett. 95, 250503 (2005^ 
[35] A. M. Childs, E. Farhi, and J. Preskill, Phys. Rev. A 65, 012322 (2002). 

[36] J. Chiaverini, D. Leibfried, T. Schaetz, M. D. Barrett, R. B. Blakestad, ,J. Britton, W. M. Itano, J. D. Jest, E. Knill, 

C. Langer, et al.. Nature 432, 602 (2004). 
[37] D. Kielpinski, V. Meyer, M.A. Rowe, C.A. Sackett, W.M. Itano, C. Monroe, and D.J. Wineland, Science 291, 1013 (2001). 
[38] Y. S. Weinstein, C. S. Hellberg, and J. Levy, Phys. Rev. A 72, 020304(R) (2005). 



